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Pathwise control of noisy bistable systems 
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Abstract 

We introduce a new method, allowing to describe slowly time-dependent Langevin 
equations through the behaviour of individual paths. This approach yields consider- 
ably more information than the computation of the probability density. The main idea 
is to show that for sufficiently small noise intensity and slow time dependence, the vast 
majority of paths remain in small space-time sets, typically in the neighbourhood of 
potential wells. The size of these sets often has a power-law dependence on the small 
parameters, with universal exponents. The overall probability of exceptional paths 
is exponentially small, with an exponent also showing power-law behaviour. The re- 
sults cover time spans up to the maximal Kramers time of the system. We apply our 
method to three phenomena characteristic for bistable systems: stochastic resonance, 
dynamical hysteresis and bifurcation delay, where it yields precise bounds on transi- 
tion probabilities, and the distribution of hysteresis areas and first-exit times. We also 
discuss the effect of coloured noise. 
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1 Introduction 

Noise is often used to model the effect of fast degrees of freedom, which are too involved to 
describe otherwise. In statistical physics and solid state physics, for instance, the influence 
of a heat bath is represented by a stochastic Glauber dynamics or a Langevin equation. 
In meteorological and climate models, the effect of fast modes (e. g. short wavelength 
modes neglected in a Galerkin approximation) is often described by noise |Ha, Ar2(| . As a 



consequence, stochastic differential equations are widely used to model systems of phys- 
ical interest, including ferromagnets |Mar|, lasers HL, HN|, neurons [Tu, glacial 



cycles | BPSV |, oceanic circulation HCef , biomolecules [ SHD |, and more. 

Among the simpler stochastic models in use is the Langevin equation with additive 
white noise 

dxt = -VV{xt,X)dt + aG{X)dWt. (1.1) 

Here F is a potential, Wt denotes a standard vector- valued Wiener process (i. e., a Brow- 
nian motion), and a measures the noise intensity. For now, we consider A as a fixed 
parameter, but below we will be concerned with situations where A varies slowly in time. 



Of course, one may be interested in situations where noise enters in a different way, for 
instance G depending on x as well, or coloured noise. 

There exist different methods to characterize the dynamics of the Langevin equation 
(|1.1| ). A popular approach is to determine the probability density p{x, t) of xt, which gives 
all information on the instantaneous state of the system. For instance, the probability 
that xt belongs to a subset T) of phase space is given by the integral of t) over T). The 
density is given by the normalized solution of the Fokker-Planck equation 

^p(x, i) = V • (Vy(x, X)p{x, t)) + V • D{X)Vp{x, t), (1.2) 

where D = {a"^ /2)GG'^ is the diffusion matrix. In particular, in the isotropic case GG^ = 



1, (1^) admits the stationary solution 



Mx) = ^e-^y(^y^\ (1.3) 

where is the normalization. For small a, the stationary distribution is sharply peaked 
around the minima of the potential. For an arbitrary initial distribution, however, the 
Fokker-Planck equation cannot be solved in general, and one has to rely on spectral 
methods, WKB approximations and the like. 

Even if we have obtained a solution of (|l.2| ), this approach still has serious shortcom- 
ings. The reason is that the probability density only gives an instantaneous picture of the 
system. If, for instance, we want to compute correlation functions such as E{xsXt} for 
< s < t, solving the Fokker-Planck equation with initial condition xq is not enough: We 
need to solve it for all initial conditions {xs,s). Quantities such as the supremum of 
over some time interval are even harder to handle. 

It is important to take into account that the stochastic differential equation ( |1.1D does 
not only induce a probability distribution of xt, but also generates a measure on the paths, 
which contains much more information. For almost any realization Wt{uj) of the Brownian 
motion, and any deterministic initial condition xq, the solution {xt{uj)}t^Q of (|1.1| ) is a 
continuous function of time (though not differentiable) . The random variable xt{uj) for 
fixed t is only one of many interesting random quantities which can be associated with the 
stochastic process. 

First-exit times form an important class of such alternative random variables, and have 
been studied in detail. If D is a (measurable) subset of phase space, the first-exit time of 
Xt from D is defined as 

TT,:=mi{t>0: xt^V}. (1.4) 

For instance, if P is a set of the form {x: V{x) ^ Vi}, containing a unique equilibrium 
point which is stable, then the distribution of r-p is asymptotically exponential, with 
expectation behaving in the small-noise limit like Kramers' time 

TKramers = e^^^i"^")/-' , where Vq := min V{x). (1.5) 

A mathematical theory allowing to estimate first-exit times for general n-dimensional 
systems (with a drift term not necessarily deriving from a potential) has been developed 



by Freidlin and Wentzell |FW]. In specific situations, more precise results are available, for 
instance the following. If D contains a unique, stable equilibrium point, subexponential 
corrections to the asymptotic expression (|1.5|) are known, even for possibly time-dependent 
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drift terms [Xz, |FJ[|. The case where contains a saddle as unique equihbrium point has 



been considered by Kifer in the seminal paper |Ki| . The situation where D contains a stable 
equilibrium in its interior and a saddle on its boundary is dealt with in [ MS1| ], employing 
the method of matched asymptotic expansions. The limiting behaviour of the distribution 
of the first-exit time from a neighbourhood of a unique stable equilibrium point as well as 



from a neighbourhood of a saddle point has been obtained by Day | Dayl| , |Day2 | 



We remark that more recently, another approach has been introduced, which mimics 
concepts from the theory of dissipative dynamical systems |CF1, |Schm , Arl| . The main 



idea is that for a given realization of the noise (i.e., in a quenched picture), paths with 
different initial conditions may converge to an attractor, which has similar properties as 
deterministic attractors. Of course, in experiments, this random attractor is only visible 
if we manage to repeat the experiment many times with the same realization of noise. 

The method of choice to study the Langevin equation ( |1.1| ) also depends on the time 
scale we are interested in. Consider for instance a one-dimensional double-well potential 
V{x) = —^ax^ + \hx'^, a,b > 0, with a barrier of height H = a?/Ab. Assume that xq 
is concentrated at the bottom x = \fajh of the right-hand potential well. If the noise 
is sufficiently weak, paths are likely to stay in the right-hand well for a long time. The 
distribution of xi will first approach a Gaussian in a time of order 

^relax = (1-6) 
C 

where c = 2a is the curvature at the bottom of the well (the variance of the Gaussian is 
approximately o"^/ (2c)). With overwhelming probability, paths will remain inside the same 
potential well, for all times significantly shorter than Kramers' time Txramers = ■ 
Only on longer time scales will the density of xi approach the bimodal stationary density 
(|l.3| ). The dynamics will thus be very different on the time scales t <C Treiax, ^reiax <C t ^ 
ramers) and t 3> ^Kramers- Random attractors can be reached only in the last regime. 
In particular, results in ||CF2 | stating that for the double- well potential V ^ the random 



attractor almost surely consists of one (random) point apply to the regime t ^> Txramers- 
In the present work, we are interested in situations where the parameter A varies slowly 
in time, that is, we consider stochastic differential equations (SDEs) of the form 

dxt = -Vy(xt, A(et)) dt + aG{\{et)) dlVj. (1.7) 

Such situations occur if the system under consideration is slowly forced, for instance by 
an external magnetic field, by climatic changes, or by a varying energy supply. Note that 
the probability density of xt still obeys a Fokker-Planck equation, but there will be no 
stationary solution in general. 

The main idea of our approach to equations of the form ( |1.7D is the following. For 
sufficiently small e and o", and on an appropriate time scale, we show that the paths {xjjt^o 
can be divided into two classes. The first class consists of those paths which remain in 
certain space-time sets, typically in the neighbourhood of potential wells (but in some 
cases, they may switch potential wells). The geometry and size of these sets depends on 
noise intensity and shape of potential. The second class consists of the remaining paths, 
which we do not try to describe in detail, but whose overall probability is small, typically 
exponentially small in some combination of e and a. In this way, we obtain in particular 
concentration results on the density of Xt without solving the Fokker-Planck equation, 
but we also obtain global information on the stochastic process, including correlations, 
first-exit times, transition probabilities and, for instance, the shape of hysteresis cycles. 
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The slow time dependence of A introduces a new time scale. If, for instance, A is 



periodic with period 1, then (1^) depends periodically on time with period Tforcing = 



Since the shape of the potential changes in time, the definitions (1.5) and ( [L.6[ ) of the 
Kramers and relaxation times no longer make sense. We can, however, define these time 
scales by 

T^Sers = e^^-/^^ i;£^ = -^, (1.8) 

where -ffmax and Cmax denote, respectively, the maximal values of barrier height and cur- 
vature of a potential well over one period. Our results typically apply to the regime 

^(min) rp „ ^(max) ,^ „n 

-'relax ^ forcing ^ Kramers' 

that is, we require that e ^ Cmax and <^ 2//inax/|loge|. The minimal curvature and 
barrier height, however, are allowed to become small, or even to vanish. 

We can describe the paths' behaviour on time intervals including many periods of the 
forcing, as long as they are significantly shorter than the maximal Kramers time. We find 
it convenient to measure time in units of the forcing period. After scaling f by a factor e, 
the relevant time scales become 

^(min) _ £ rj. ^(max) _ 2H,^,^/a^ ^^ ^n\ 

-'relax ~„ ' -'forcing — -'Kramers ■ \^-^^) 

t-max 

When scaling the Brownian motion, we should keep in mind its diffusive nature, which 
implies that in the new units, its standard deviation grows like \ftfe. Equation (^]^) thus 
becomes 

Axt = --VV{xt, X{t)) dt + ^ G{X{t)) dWt. (1.11) 
£ ye 

In the deterministic case u = 0, we will sometimes write this equation in the form ex = 
— Vy(3;, A(t)), which is customary in singular perturbation theory. Throughout the paper, 
we will assume that V, G and A are smooth functions, that V{x, A) grows at least like 
for large ||x||, and that the matrix elements of G, as well as A, have uniformly bounded 
absolute values. A large part of the paper is devoted to one-dimensional systems, and we 
come back to the multidimensional case in Section ^. 

The aim of this paper is to illustrate our methods by applying them to a number of 
physically interesting examples. We thus emphasize the conceptual aspects, and refrain 



from giving mathematical details of the proofs, which will be presented elsewhere [ BGl 
BG2|, pG| 



We start, in Section by explaining the basic ideas in the simplest situation, where 
most paths remain concentrated near the bottom of a potential well. We also briefly 
describe the dynamics near a saddle. The three subsequent sections are devoted to more 
interesting cases in which paths may jump between the wells of a double-well potential. 
Note that we restrict our attention to double-well potentials only to keep the presentation 
simple, but the dynamics in more complicated multi-well potentials can be described by 
the same approach. 

In Section ^, we discuss the phenomenon of stochastic resonance, where noise allows 
transitions between potential wells which would be impossible in the deterministic case. 
We compute the threshold noise intensity needed for transitions to be likely, and show 
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Figure 1 . Example of a periodically forced double- well potential. In the limit of infinitely 
slow forcing, an overdamped particle tracks the bottom of potential wells. However, its 
position is not determined by the instantaneous value of the forcing alone. The system 
describes a hysteresis cycle, whose shape depends on the frequency of the forcing. Noise 
may kick the particle over the potential barrier, and thus influences the shape of the 
hysteresis cycle. 

that most paths are close, in a natural geometrical sense, to a periodic function. This 
provides an alternative quantitative measure of the signal's periodicity to the commonly 
used signal-to-noise ratio. 

Section ^ is devoted to hysteresis, which is also characteristic for forced bistable sys- 
tems. For sufficiently strong forcing, one of the two potential wells disappears, causing 
trajectories to switch between potential wells even when no noise is present (Figure ||). 
As a result, even in the adiabatic limit, the instantaneous value of the parameter A does 
not suffice to determine the state of the system: Solutions follow hysteresis cycles. In 
absence of noise, their area is known to scale in a nontrivial way with the frequency of the 
forcing. Our methods allow us to characterize the distribution of the random hysteresis 
area for positive noise intensity. In particular, we show that for noise intensities above an 
amplitude-depending threshold, the typical area no longer depends, to leading order, on 
amplitude or frequency of the forcing. 

In Section ^, we consider the effect of additive noise on systems with spontaneous 
symmetry breaking, i.e., when the potential transforms from single to double well. In 
the deterministic case, solutions are known to track the saddle for a considerable time 
before falling into one of the potential wells, a phenomenon known as bifurcation delay. 
We characterize the effect of additive noise on this delay, and on the probability to choose 
one or the other potential well after the bifurcation. Furthermore, we give results on the 
concentration of paths near potential wells which allow, in particular, to determine the 
optimal relation between speed of parameter drift and noise intensity for an experimental 
determination of the bifurcation diagram. 

Finally, Section ^ contains some generalizations. We first discuss analogous results to 
those of Section ^ for multidimensional potentials. This formalism allows us to treat the 
effect of the simplest kind of coloured noise, given by an Ornstein-Uhlenbeck process, in a 
natural way. We conclude by discussing the dependence of previously discussed phenomena 
on noise colour. 

Acknowledgements: 

We are grateful to Anton Bovier for helpful comments on a preliminary version of the 
manuscript. N. B. thanks the WIAS for kind hospitality. 
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2 Near wells and saddles 



We start by discussing situations in which the noise intensity is sufficiently small, compared 
to the depth of a given potential well, for paths to remain concentrated near the bottom 



of the well during a long time interval. In Section 2.1, the solvable linear case is used to 
compare the information provided by the Fokker-Planck equation and by the pathwise 
approach. In Section |2.2|, we show that our method naturally extends to the nonlinear 



case. We briefly describe the dynamics near a saddle in Section 2.3 



2.1 Linear case 

It is instructive to consider first the case of a linear force (that is, of a quadratic potential), 
which can be solved completely. A general one-dimensional, time-dependent quadratic 
potential can be written as 

y(x,^) = ^c(^)(x-x*(t))^ (2.1) 

where x*(t) is the location of the potential minimum, and c(t) is the curvature of the 
potential. The SDE ( [1.7| ) takes the form 

dxt = -a%t) {xt - x*{t)) dt + -^g{t) dWt, (2.2) 
e ye 

where a*{t) = —c{t). Throughout this paper, we will use a* to denote the linearization of 
the force at an equilibrium point x* , with a* < if the equilibrium is stable, and a* > 
if it is unstable. In this subsection and the following one, we consider the stable case, and 
assume that a*{t) ^ — oq for all t under consideration, where oq is a positive constant. 
For simplicity, we shall assume that the functions x*{t) and a*{t), as well as g{t), are 
real-analytic. 

Let us first investigate the probability density p{x,t) of Xf. The Fokker-Planck equa- 
tion being linear, it can be easily solved. Assume for simplicity that the distribution of xq 
is Gaussian, with expectation E{xo} and (possibly zero) variance Var{xo}. (In addition, 
we always assume the initial distribution to be independent of the Brownian motion.) 
Then xt has a Gaussian distribution for any t > 0, with density 

p(x,t) = ^ expl-^ \-^}, 2.3) 

^ ^ V2^Var{xa "^l 2YaT{xt} J' ^ ^ 

where expectation and variance of xt obey the ODEs 

^E{xt} = ^a*{t) {E{xt} - x*{t)) (2.4) 

Var{xt} = -a*{t) Var{xt} + —git)^. (2.5) 
dt e e 

Note that Ejxf} coincides with the deterministic solution xf^* of Equation (|2.2| ) for a = 0, 
with initial condition Xq'^^ = K{xo}: 

E{xt} = * = 4"* e"*W/" — / e"*(*'^)/^ a*(s)x*(s) ds, (2.6) 

£ Jo 
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where we use the notations 

a*{t,s)= I a*{u)du, a*{t) = a*{t,0). (2.7) 

J s 

Our stabUity assumption a*{t) ^ — oq Vt imphes that a*{t,s) ^ —ao{t — s) for t > s. 
Hence the first term on the right-hand side of ( |2.6| ) decreases exponentially fast: It is at 
most of order e after time e|loge|/ao, at most of order after time 2e|loge|/ao, and so 
on. 

We expect xf'^^ to follow adiabatically the slowly drifting bottom of the potential well. 
To make this apparent, we evaluate the second term on the right-hand side of ( |2.6| ) by 
integration by parts: 

_1 r QO.*(t,s)/e a*(s)x*{s) ds = x*{t) - x'^(O) e"*W/^ - T e"*(*'^)/^ x*{s) ds. (2.8) 
£ Jo Jo 



By successive integrations by parts, we find that the general solution of ( |2.4| ) can be 
written as 

E{xt} = xf * = xf' + (x[J"* - x^"*) e"*W/S (2.9) 
where xf^^ is a particular solution of ( |2.4D , admitting the asymptotic expansion]^ 
.det _ , , .2 1 d fx^it) 



xf^ = x*{t) + e—^ + e'—-^{—fi]+--- (2.10) 

Since a*(t) is negative, x^'^* tracks the bottom x*(t) of the potential well with a small lag: 
^det ^ x*(t) if x*(t) moves to the right, and x^*^* > x*(t) if x*(t) moves to the left. The 
particular solution ( |2.1[1| ) is called adiabatic solution or slow solution. Relation (2.9) ex- 
presses the fact that all solutions of (|2.4| ) are attracted exponentially fast by the adiabatic 
solution xf*^*. 

The variance of xt can be computed in a similar way. The solution of ( |2.5| ) is given by 

Var{xa = Var{xo}e2°*W/^+— / e^"*^*'")/" ^(s)^ ds. (2.11) 

£ JO 

The behaviour of the variance is very similar to the behaviour of the deterministic solution 



2.f:): The initial condition Varjxo} is forgotten exponentially fast, and in analogy with 
2.E ) and ( ^.10 ), we can write 



Var{xt} = v{t) + (Var{xo} - v(0)) e^"*^/^ (2.12) 
where v{t) is a particular solution of (p.5|), which admits the asymptotic expansion 



vit) 



(2.13) 



2|a*(t)| 

The fact that xt has a Gaussian distribution implies in particular that for any t > 0, 



\xt-xf''\ 



= 2 / , p(xf * + y, t) dy ^ e-^'l^ . (2.14) 

Jhy/vMyt} 



^The asymptotic series does not converge in general, but it admits expansions to any order in e, with 
a remainder which can be controlled. 
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Hence, the distribution of Xt is concentrated in an interval of width y^2 Varjyt} around 
Xi^^, which behaves asymptotically like crg{t)/ \/\a*{t)\ by ( |2.12 ) and ( ^.13 ). In words, 



the spreading of xt is proportional to the noise intensity and inversely proportional to 
the square root of the curvature of the potential: Flatter potentials give rise to a larger 
spreading of the distribution. 

Up to now, we have only studied the probability density of Xf. However, even if the 
density is concentrated near the bottom of the well at all times, this does not exclude that 
the path {a^sjo^s^^ makes occasional excursions away from x*. From now on, we consider 
the initial condition xq as deterministic, so that the path depends only on the realization 
of the Brownian motion. The solution of the SDE (12.21) can be written as 



Xt 



f'' + yu yt = ^ [\'^*^'''^/'g{s)dWs. (2.15) 

Jo 



The process yt is a generalization of an Ornstein-Uhlenbeck process, with time-dependent 
damping and diffusion. Ideally, we would like to estimate the probability that the path 
leaves a strip of (time-dependent) width proportional to y^Varjxt}, and centred at x^®*. 
This turns out to be difficult because the variance may change quickly for t very close 
to 0, due to the first term on the right-hand side of ( p.ll]) . To avoid these technical 
complications, we use a strip of width proportional to s/^it), defined by 



Bih) = {{x,t): \x-xr^\ <h^/^}. (2.16) 

Note that B{h) coincides, up to order e, with the set of points where V{x, t) is smaller 
than V{xf'^^,t) + {^hag{t))'^ . Showing that the path {x^lo^s^t is likely to remain in B{h) 
is equivalent to showing that the first-exit time of Xs from B{h), defined by 

rB(/,) =inf{O0: {xs, s) ^ B{h)} , (2.17) 

is unlikely to be smaller than t. In fact, the following probabilities are equivalent (the 
superscripts refer to the initial condition): 



<t} =P°'""{35G [0,i): {xs,s)^B{h)} 

( U _ ^det| 



sup ^ ' ^ (2.18) 



The following result is a straightforward consequence of standard exponential bounds on 
the supremum of stochastic integrals, extended to integrals as appearing in ( |2.15| ). The 
proof, given in [BGl, Proposition 3.4] for constant g, also applies here.0 



Proposition 2.1. For all t and h > 0, 



IP°'"M^BW<t}^C(t,e)e-'^'^ , (2.19) 

where 



\n*(f)\ 1 
C(t,e) = J — ^ + 2 and K=--0{e). (2.20) 



^The generalization of the proof to g bounded away from zero is trivial, but the result also holds if g 
vanishes. In fact, a sufficient condition is that v{t)/v{s) = 1 + 0(e) whenever t — s = 0(e^), which can be 



checked using the asymptotic expansion (2.13) 
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Figure 2. A sample path xt of the Unear equation ( |2.2D for a*{t) — —A + 2siii(47rt), 
x*{t) = sin(27rt) and g{t) = 1. Parameter values are e = 0.04 and a = 0.025. The shaded 
region is the set B{h) for h — 3, centred at the deterministic solution xf'^^ starting at the 



same point as Xt- The width of B{h) is of order ha/ ^a*{t). Proposition 2.1 states that 
the probability of Xt leaving B{h) before times of order 1 decays roughly like e"'* 

Note that only the prefactor C{t, e) is time-dependent. The exponential factor e"'^'*^ is 
small as soon as /i ^ 1, so that the paths {x^jo^s^t are concentrated in a neighbourhood of 
order ^/v{s) of the deterministic solution up to time t (see Figure |2|) . More precisely, paths 
are unlikely to leave the strip B{h) of width h^/^ before time t provided h? » logC(t,e). 

The bound (|2.19| ) is useful for times significantly shorter than Kramers' time, which is 
of order e^^l'^ (to reach points where the potential has value i?). For longer times, the 
prefactor C(t, e) becomes sufficiently large to counteract the term e"'^'* for any reasonable 
/i, which is natural, as we cannot expect paths to remain concentrated near x^''* on such 
long time scales. On polynomial time scales of order (y~^ ^ however, large excursions are 
very unlikely. 



The estimate (2.19) has been designed to yield an optimal exponent for noise intensities 
scaling like a power of e. We do not expect the prefactor C{t,e) to be optimal, but for times 
and noise intensities polynomial in e, it leads to subexponential corrections. However, if 



we do not care for the precise exponent, (2.20) can be replaced by 



C(t,e) = ^^^ + 2 and k>0. (2.21) 

The denominator e in C is due to the fact that we work in slow time. 

2.2 Nonlinear case 

We consider now the motion in a more general, nonlinear potential of the form 

V{x,t) = ic(t)(x-x*(t))' + 0((rc-x^(t))'), (2.22) 

admitting a local minimum at x*{t). As before, we assume that the curvature c{t) is 
bounded below by a positive constant for all times t. We do not exclude, however, that V 
has other potential wells than the one at x*(t). 

The probability density can no longer be computed exactly in general, although it seems 
plausible that a distribution initially concentrated near x*{0) will remain concentrated near 
x*{t), on a certain time scale. In fact. Proposition |2.1| naturally extends to the nonlinear 
case. 

Consider first the deterministic case cr = 0. A result due to Gradstem and Tihonov 
|Gr, Ti], which is related to the adiabatic theorem of quantum mechanics, states that 
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• there exists a particular solution xf^^ of the deterministic equation ex = —dxV{x, t) 
tracking the bottom of the potential well at a distance of order e; 

• any solution xf*^* starting in a neighbourhood of x*(0) (in fact, inside the potential 
well) approaches x^^^ exponentially fast in t/e. 

Let us fix a deterministic initial condition xq = Xq^* such that xf^* is attracted by x^^^. 
We introduce the notations 



a{t) 



dx"^ 



(xf^'^t), a{t,s)= I a{u)du and a{t)=a{t,0) ( 

J s 



2.23) 



for the curvature of the potential at xf"^^ and the analogue quantities to (|2.7| ). Note that 
Tihonov's result implies that a{t) asymptotically approaches — c(i) + 0(e). The difference 
yt = xt — xf^^ satisfies an SDE of the form 



1 



dyt = - [a{t)yt + t)] dt + ^g{t) dWt, 

e ye 



(2.24) 



where b{y,t) = 0{y'^) describes the effect of nonlinearity. We define again a strip B{h) as 
in (|2l^ ), with 



v{t) = v{f)) e2"WA + 



a 



9{s) 



ds. 



(2.25) 



Our results work for any v{0) larger than a positive constant independent of e. A conve- 
nient choice is v{0) = (T^(7(0)^/(2|a*(0)|): Then the fact that xf'^^ approaches exponentially 
fast a neighbourhood of order e of x*{t) implies that 



v{t) 



a 



9{t? 



2\a*{t)\ 



+ 0{e) + 0(1x0 - x*{{))\ e-'^'"^** */^) 



(2.26) 



Proposition 2.1 generalizes to 



Theorem 2.2. There exists a constant ho, independent of a and e, such that for all 
h ^ ho/a, 



where 



C{t,e) 



\a{t)\ 



l)0,a;'o 



+ 2 



{rB(h) < t} ^ C(t,e)e 



and 



_ _ o{e) - 0{ah). 



(2.27) 



(2.28) 



The interpretation is the same as in the linear case: Paths are concentrated, for times 
significantly shorter than Kramers' time, in a strip of width proportional to \/v{t) around 



X 



det 



The proof is identical to the one of [ BGl , Theorem 2.4]. The main idea is to show 
that if the solution of the equation linearized around x^^* remains in a strip B{h), then 
the solution of the nonlinear equation ( 2.24 ) almost surely remains in the slightly larger 
strip B{h[l + 0{ah)]). 

The main difference between the nonlinear and the linear case is the condition h ^ 
ho/a, which stems from the requirement that the linear term a{t)yt in Equation ( p. 24 ) 
should dominate the nonlinear term b(yt,t) for all (xf , t) € B{h). Because of this condition, 
the result ( 2.27] ) is useful for cr^ <C k/iq/ log C(t, e). It is, however, possible to derive bounds 
for larger deviations under additional assumptions on the potential: 
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Proposition 2.3. Assume that there are constants Lq > 0, K > and n ^ 2 such that 

dV 

x—{x,t)-^K\x\'' (2.29) 
ox 

whenever \x\ ^ Lq and t ^ 0. Then there exist constants C, k > such that 

,o,xoj g^^p i^^i ^l] C(- + l) e-'^^"/'"' (2.30) 



for all t ^ 0, L ^ Lq and \xo\ ^ Lq/2. 



This result is a generalization of | BG3 , Proposition 4.3], where the case n = 4 was 
treated. 



We remark in passing that the bounds ( ^.27 ) and ( 2.29 ) are sufficient to provide es 



timates on the moments of the distribution of xt, without solving the Fokker-Planck 
equation (c. f. [ BG3| , Corollary 4.6]): 



Corollary 2.4. Assume that ( 2.29| ) holds for some n ^ 2. Then 



E°'^»{|xt - xf^Y''} ^ (2fc - iy.lM''v{t)'' (2.31) 

for some constant M, all integers k and all t ^ 0, provided a ^ Co/log(l + t/e) for a 
sufficiently small constant cq. 

Note that the Cauchy-Schwarz inequality immediately implies bounds on odd moments 
as well. The bounds on the moments are those of a Gaussian distribution with variance 
Mv{t), even if the potential V has multiple wells. The reason is that on the time scale 
under consideration, solutions of the SDE do not have enough time to cross a potential 
barrier and reach another potential well. 

2.3 Escape from a saddle 

Assume now that the potential V{x, t) admits a saddle at x*{t) for all times under consid- 
eration. In the deterministic case, a particular solution xf'^^ is known to track the saddle 
at a distance of order e, separating the basins of attraction of two neighbouring potential 
wells. Trajectories starting near x^*^* will depart from it exponentially fast, but if the 
initial separation | xq — Xq'^^ \ is exponentially small, the time required to reach a distance 
of order one from the saddle may be quite long. 

Noise will help kicking xt away from x^^^, and thus reduce the time necessary to leave 
a neighbourhood of the saddle. In order to describe this effect, we consider the deviation 
yt = Xt — , which satisfies the SDE 

dyt = - [a{t)yt + b{yt, t)] dt + -^g{t) dWt, (2.32) 
£ ye 

where a{t) ^ > is the curvature of the potential at x^'^^, and b{y,t) = 0{y'^). We 
now describe the dynamics in a small neighbourhood of the adiabatic solution tracking 
the saddle, where diffusion prevails over drift^, defined by 

B{h) = [{x,t):\x-xr\<^^^\ (2.33) 



possible way to compare the influence of drift and noise is by Ito's formula, which yields, for 6 = 0, 
A{yt) = (l/e)[2a(t)i/j + (j'^g{t)^] dt + {2a / y/e)g{t)yt dWt- The expression in brackets is dominated by the 
noise term a^g{t)^ for all {xt,t) in B{1), while the deterministic term 2a{t)y^ prevails outside this domain. 



11 



The following result, which is proved in the same way as |BG1, Proposition 3.10], shows 
that the first-exit time 'rg(/i) of xt from B{h) is likely to be small. 



Theorem 2.5. Assume that g is bounded below by Le for a sufficiently large constant L. 
Then for all h ^ 1 and all initial conditions (2;o,0) G B{h), 

>t}^ Cexp{-~l^a{s)ds], (2.34) 

where C, k are positive constants. 

This result shows that xt will leave a neighbourhood of size ag{t)/ y'2a{t) of xf'^^ 
typically after a time of order e/a(0). Once this neighbourhood has been left, the drift 
term starts prevailing over the diffusion term, and one can show, (although this is not 
trivial,) that the typical time needed to leave a neighbourhood of order one of the saddle 
is of order e|log a\. We will state a similar result in Section ^ when discussing the dynamics 
after passing through a pitchfork bifurcation point. 



3 Stochastic resonance 



Up to now, we have considered situations in which the potential has bounded curvature 
near its (isolated) extrema, so that for sufficiently small noise intensities and not too long 
time scales, most paths are concentrated near the bottom of the well they started in. 

Not surprisingly, interesting phenomena occur when the condition on the curvature is 
violated. Two cases can be considered: 

• Avoided bifurcation: The potential well becomes flatter, but the curvature does not 
vanish completely; sufficiently strong noise, however, may drive solutions to another 
potential well. This mechanism is responsible in particular for the phenomenon of 
stochastic resonance. 

• Bifurcation: The curvature at the bottom of the well vanishes, say at time 0; for 
t > 0, new potential wells may be created (e. g. pitchfork bifurcation) or not (e. g. 
saddle-node bifurcation). 

We will discuss the possible phenomena in the case of a Ginzburg-Landau potential 

V{x,t) = ^x^-^^,{t)x^ -X{t)x. (3.1) 

However, the precise form of the potential and the fact that it has at most two wells are 
not essential. 

The potential (3.1) has two wells if 27A^ < and one well if 27A^ > Crossing 



the lines 27A^ = 4/i^, jj, > 0, corresponds to a saddle-node bifurcation, and crossing the 
point A = ;U = corresponds to a pitchfork bifurcation. Equilibrium points are solutions 
of the equation x^ — fi{t)x — X{t) = 0; we will denote stable equilibria by x^-, and the 
saddle, when present, by Xq. 

In this section, we investigate situations with avoided bifurcations, in which the po- 
tential always has two wells, but the barrier between them becomes low periodically. 
Bifurcation phenomena will be discussed in the next two sections. 
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Figure 3. A sample path of the equation of motion (3.2) in an asymmetrically perturbed 
double-well potential. Parameter values are e = 0.0025, oq = 0.005 and a = 0.065, which 
is just above the threshold for transitions to be likely. The upper and lower full curves 
show the locations of the potential wells, while the broken curve marks the location of the 
saddle. 



3.1 The mechanism of stochastic resonance 

Let us consider the case where ^ is a positive constant, say /u = 1, and A(t) varies period- 
ically, say A(t) = — Acos(27rt). If |A| < Ac = 2/(3-v/3), then F is a double-well potential. 
We thus assume that j4 < Ac- 

In the absence of noise, the existence of the potential barrier prevents the solutions 
from switching between potential wells. If noise is present, but there is no periodic driving 
i^A = 0), solutions will cross the potential barrier at random times, whose expectation is 
given by Kramers' time ee^^/'^ , where H is the height of the barrier {H = 1/4 in this 
case) . 

Interesting things happen when both noise and periodic driving A(t) are present. Then 
the potential barrier will still be crossed at random times, but with a higher probability 
near the instants of minimal barrier height (i.e., when t is integer or half-integer). This 
phenomenon produces peaks in the power spectrum of the signal, hence the name stochastic 
resonance (SR). 

If the noise intensity is sufficiently large compared to the minimal barrier height, 
transitions become likely twice per period (back and forth), so that the signal xt is close, 
in some sense, to a periodic function (Figure |3|). The amplitude of this oscillation may be 
considerably larger than the amplitude of the forcing A(t), so that the mechanism can be 
used to amplify weak periodic signals. This phenomenon is also known as noise-induced 



synchronization |SNA| , |NSAS| . Of course, too large noise intensities will spoil the quality 
of the signal. 

The mechanism of stochastic resonance was originally introduced as a possible ex- 
planation of the close-to-periodic appearance of the major Ice Ages [ BSV , BPSV| . Here 
the (quasi-)periodic forcing is caused by variations in the Earth's orbital parameters (Mi- 
lankovich factors), and the additive noise models the fast unpredictable fluctuations caused 
by the weather. Meanwhile, SR has been detected in a large number of systems (see for 
instance [MW, WM, pHM ] for reviews), including ring lasers, electronic devices, and even 
the sensory system of crayfish and paddlefish | Nfc |. 

Despite the many applications of SR, its mathematical description has remained in- 
complete for two decades, although several limiting cases have been studied in detail. The 
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first approaches considered either potentials that are piecewise constant in time jBSV], 
or two-level systems with discrete space [ [ET] , pVIcNW ], Continuous time equations have 
been mainly investigated through the Fokker-Planck equation, using methods from spec- 
tral theory | Fox , pHl | or linear response theory |JH|. The main contribution of these 
approaches is an estimation of the signal-to-noise ratio (SNR) of the power spectrum, as 
a function of the noise intensity. The SNR is one of the possible quantitative measures 
of the signal's periodicity, and behaves roughly like /cj^, which is maximal for 

= H/2. In [ |MS2 |, an action functional is used to extend Kramers' result to the case of 
small-amplitude forcing. 

A description of individual paths has been given for the first time in Freidlin's recent 
paper |Frl]. His results apply to a general class of n-dimensional potentials, in the case 
where the period 1/e of the driving scales like Kramers' time ^^1" . The fact that the 
minimal barrier height K is considered as constant while a tends to zero, implies that the 
results only hold for exponentially long driving periods. As quantitative measure of the 
signal's periodicity, the L^'-distance^ between paths {xtjf^o and a deterministic, periodic 
limit function is used. This limit function simply tracks the bottom of a potential 
well, and jumps to the deeper well each time the potential barrier becomes lowest. The U*- 
distance is shown to converge to zero in probability as a goes to zero. However, Freidlin's 
techniques do not yield estimates on the speed of this convergence, or its dependence on -p. 

Our techniques allow us to provide such estimates for one-dimensional potentials, with 
a more natural distance than the L^-distance: In fact, we simply use a geometrical distance 
between paths and the limit function, considered as curves in the (t, x)-plane. The analysis 
given below also includes situations in which the minimal barrier height becomes small in 
the small- noise limit. 



3.2 Pathwise description 

For the Ginzburg-Landau potential ( |3.1] ), = 1 and \(t) = — ^cos(27rt), the SDE takes 
the form 

dxj = — \xt — xi — A cos {2^t)\dt + —<iWt. (3.2) 

£ ye 

We assume that ^ < Ac, so that there are always two stable equilibria at and a 

saddle at XQ{t). We introduce a parameter ao = Xc — A which measures the minimal 
barrier height: At t = 0, the barrier height is of order Og^^ for small ao, and the distance 
between x'^ and the saddle at Xq is of order ^/oq. At t = ^, the left-hand potential well 
at X* is likewise close to the saddle. In order for transitions to become possible on a time 
scale which is not exponentially large, we allow oq to become small with e. 

Assume that we start at time = — 1/4 in the basin of attraction of the right-hand 
potential well. Results from Section ^ show that transitions are unlikely for t <^ 0. Also, 
for <^ t <ti 1/2, paths will be concentrated either near x^ or near . This allows us to 
define the transition probability as 

i^trans = P*°'"°{xt, < O}, to = "1/4, h = 1/4- (3.3) 

The properties of Ptrans do not depend sensitively of the choices of to and ti, as long as 
— 1/2 ^ to <C <C ti <C 1/2. Also the level can be replaced by any level lying between 

''The L''-distance between xt and (f){t) is the integral of \\xt — ||^ over a given time interval (raised 
to the power 1/p). Note that in contrast to a small supremum norm, a small L'^-norm of xt — (t>{t) does not 
exclude that Xt makes large excursions away from ^(t), as long as these excursions are sufficiently short. 
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Figure 4. Sample paths of Equation (3.2) in a neighbourhood of time 0, when the barrier 
height is minimal, for two different noise intensities. Full curves mark the location x\{t) 
of the right-hand potential well, broken curves the location XQ{t) of the saddle. Parameter 
values are e = 0.0125, ao = 0.002 and (a) a = 0.012, (b) a = 0.07. In case (a), the 
path remains in the set B{h), shown here for /i = 3, which is centred at the deterministic 
solution a;^°'. In case (b), Xt remains in B{h) as long as the width of B{h) is smaller than 
the distance between a;^''* and the saddle, that is, for t ^ — cr^/'^. The path jumps to the 
left-hand potential well during the time interval 



x*_{t) and x'X{t) for all t. Denoting by a V 6 the maximum of two real numbers a and b, 
our main result can be formulated as follows. 

Theorem 3.1 ([^G|, Theorems 2.6 and 2.7]). For the noise intensity, there is a thresh- 
old level (Tc = (ao V e)^/^ with the following properties: 

1. If (J < CTc, then 

Ptr.ns^^e"''"''"" (3.4) 

— 1/3 

for some C,k > 0. Paths are concentrated in a strip of width a/{y'\t\ V ac ) around 
the deterministic solution tracking x*^{t) (Figure^a). 

2. If a > ac, then 

Ptrans>l-Ce-'^^'''/(^|l°S'^l) (3.5) 

for some C, k > 0. Transitions are concentrated in the interval [— o"^/^, cr^/^]. More- 
over, for t ^ —a^l^, paths are concentrated in a strip of width a / \/\t\ around the 
deterministic solution tracking x*i^{t), while for t ^ a"^^^ , they are concentrated in a 
strip of width cr/\/t around a deterministic solution tracking x* (t) (Figure^). 

The crossover is quite sharp: For o" <C cTc, transitions between potential wells are very 
unlikely, while for a ^ ac, they are very likely. By "concentrated in a strip of width 
w" , we mean that the probability that a path leaves a strip of width hw decreases like 
g-K/i £qj. gQjj^g K > 0. The "typical width w" is our measure of the deviation from the 
deterministic periodic function, which tracks one potential well in the small-noise case, 
and switches back and forth between the wells in the large-noise case. 



Theorem 3.1 implies in particular that for the periodic signal's amplification to be 



optimal, the noise intensity a should exceed the threshold ac- Larger noise intensities will 
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increase both spreading of paths (especiahy just before they cross the potential barrier) 
and size of transition window, and thus spoil the output's periodicity. 

It has been proposed to identify stochastic resonance and synchronization with a phase- 
locking mechanism [ |NSAS |, but the definition of a phase remained problematic. The fact 



that the majority of paths are contained in small spacer-time sets makes it possible to 
associate a random phase with those paths. For instance, when a > ctc, most paths circle 
the origin of the (A, x)-plane counterclockwise (see Figure ^c). For all paths not containing 
the origin, one can define a random phase ^pt by tan((/9t) = xt/X{t). Although (ft fluctuates, 
it is likely to increase at an average rate of 27r per cycle. 

Part of our results should appear quite natural. If oq > £, the threshold noise level 
(Tc = Hq^^ behaves like the square root of the minimal barrier height, which is consistent 
with the SNR being optimal for o"^ = H/2. However, ac saturates at e^/^ for all oq ^ 
e. Hence, even driving amplitudes arbitrarily close to Ac cannot increase the transition 
probability. This is a rather subtle dynamical effect, mainly due to the fact that even if 
the barrier vanishes at i = 0, it is lower than e^/^ during too short a time interval for 
paths to take advantage. The situation is the same as if there were an "effective potential 
barrier" of height proportional to a^. 

Another remarkable fact is that for a > ac, neither the transition probability nor 
the width of the transition windows depend on the driving amplitude to leading order. 
In the remainder of this subsection, we are going to explain some ideas of the proof of 



Theorem 3.1, which will hopefully clarify some of the above surprising properties. 

The first step is to understand the behaviour of the solutions of ( |3.2D in the determin- 
istic case a = 0. This problem belongs to the field of dynamical bifurcations, the theory 
of which is relatively well developed. We follow here a framework allowing to determine 



scaling laws of solutions near bifurcation points, which has been presented in [BK|. The 
main idea is that when approaching a bifurcation point {t*,x*), the distance between the 
adiabatic solution and the static equilibrium branch scales like e/\t — t*\P for t ^ — e'^ , 
and like for |t — t*| ^ e'^ . The rational numbers p, v and \x = \ — pv are universal 
exponents, which can be deduced from the Newton polygon of the bifurcation point. 

Tihonov's theorem implies that away from the avoided bifurcation point at t = 0, 
solutions x\^^ of the deterministic equation track the equilibrium branch x\{t) adiabat- 
ically at a distance of order e. It is thus sufficient to understand what happens in a 
small neighbourhood of the almost-bifurcation. Note that if A = Ac, the right-hand well 
and the saddle merge at x = l/\/3- It is thus helpful to consider the translated variable 
ydet _ ^det _ i/y^^ which obcys the differential equation 

eij = ct^ — \/3y^ + ao + higher order terms, c = 27r^Ac. (3-6) 



Consider the "worst" case ag = 0. Then the right-hand side of ( |3.6D describes a trans- 
critical bifurcation between the equilibria q ~ ±y^3^"'^/^|t| + ©(t'^). One shows that 
the solution yf*" tracks y\{t) at a distance scaling like for |t| ^ y^, and like for 
\t\ ^ \fe. In fact, yf^^ never approaches the saddle closer than a distance of order -y/e, 
which is because the term — \/3y^ only dominates during a short time interval of order 

The same qualitative behaviour holds for < oq ^ e. For ao > e, one can show that 
x^^"^ tracks x*^(i) at a distance never exceeding Oiej ^fa^. Since a;^(0) — x* (0) is of order 
y/oo, x^'^* never approaches the saddle closer than a distance of order 

Let us now consider the random motion near x^'^^ for o" > 0. We denote, as usual, 
by a(t) the curvature of the potential at the deterministic solution x^'^*. By (|3.6|), a{t) 
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behaves, near t = 0, hke — yf*^*, which we know to behave hke —{\t\ V for oq > £, and 



Uke —{\t\ V ^/e) for ao ^ S- It turns out that Theorem 2^ can be extended to the present 
situation, to show that paths are concentrated in a strip around x^^*. The width of this 
strip is again related to the standard deviation of a linearized process, and behaves like 

a a 



(3.7) 

However, this property only holds under the condition that the spreading is always smaller 
than the distance between x^^^ and the saddle, which scales like |f| V ctc^^ ■ We thus have 
to require a < V CTc, so that 

• if cr < (7c, the condition is always satisfied, and thus ( |3.4D follows from the generaliza- 



tion of Theorem 2.2 with h = ac/cr (Figure ^a); 
• if o" > (7c, the condition is only satisfied for t ^ —a"^/^ (Figure |^). 

It thus remains to understand what happens for t ^ —a^l"^ \i a > a^- Here the main 
idea is that during the time interval [—(7^/^, (7^/^], the process xt has a certain number of 
trials to reach the saddle. If xt reaches the saddle, it has roughly probability 1/2 to move 
in each direction. If it moves far enough towards the left, it will most probably fall into 
the left-hand potential well, and is unlikely to come back for the remaining half-period. If 
it moves to the right, it has failed to overcome the barrier, but can try again during the 
next excursion. 

One can define a typical time At for an excursion as the time needed for a path 
starting near to reach and overcome the barrier with non-negligible probability, say 
with probability 1/3. One can show, by comparison with suitable linearized processes, 
that this typical time is determined by the condition 

|a(s) I ds = const e|log(7|. (3.8) 

To obtain this, one first checks that the curvature of the potential is the same, up to sign 
reversal, at the adiabatic solutions tracking the bottom of the well and the saddle. To 



overcome the saddle, the integral in (|3.8| ) must be of order e, similarly as in Theorem 2.5 . 
The factor |log(7| is needed to reach a distance of order 1 from the saddle. Now the 
maximal number N of trials during the interval [— (7^/'^, cr^/^] is given by 

U((j2/3 _(J2/3^| ^4/3 

N = const — — ; ^ const — r. (3.9) 

e|log(7| e|log(7| 

Finally, the Markov property implies that the probability not to overcome the saddle 
during N trials is bounded by 

-) =exp|-iVlog(-)|^exp|-const^|, (3.10) 

which proves (|3.5| ). An important point to note is that the transition probability is not 
determined by the curvature of the potential at the saddle, but by the curvature at the 
deterministic solution tracking the saddle, which may be larger for small oq. 
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Figure 5. A sample path of the equation of motion (3.11) in a periodically modulated 



symmetric double-well potential. Parameter values are e = 0.005, aq = 0.005 and a ~ 
0.075, which is above the threshold for transitions to occur with probability close to 1/2. 
The upper and lower full curves show the locations of the potential wells, while the broken 
line marks the location of the saddle. 

3.3 Symmetric potentials 

Another case of interest is the Ginzburg-Landau potential ( ^.l] ) with A = and = ao + 

1 — cos(27rt), ao > 0. Then V[x^t) is always symmetric, with minima 

and a barrier height becoming small at integer times. The associated SDE is 

dxt = - [(ao + 1 - cos{2-Kt))xt - xf\ dt + dWt. (3.11) 
£ ye 

As before, transitions between potential wells are most likely when the barrier is lowest. 
We can thus define a transition probability as in ( |3.3D , with — l<Cto^O<Cti<Cl. 
We again assume that the process starts in the right-hand well. In this case, the result 
corresponding to Theorem 3.1 is 



Theorem 3.2 ([ pG^ , Theorems 2.2-2.4]). Th ere is a threshold noise level Uc = oq V e > 



with the following properties: 
1. If (J < cTq, then 



Ptr.n.^je--'^'/'^' (3.12) 



for some C,n > 0. Paths are concentrated in a strip of width o'/{\t\ V ^/^^) around the 
deterministic solution tracking x^{t). 
2. If a > ac, then 

Ptrans^^-Ce-«'^''V(.|'°S-I) (3.13) 

for some C, k > 0. Transitions are concentrated in the interval [—y/o', V^]- Moreover, 
fort ^ —\f^, paths are concentrated in a strip of width (j/\t\ around the deterministic 
solution tracking x* (t), while for t ^ ^/a, they are concentrated in a strip of width 
a/t around a deterministic solution tracking either x^{t) or x*_{t). 

The main difference with respect to the previous case is that due to the symmetry, Ptrans 
can never exceed 1/2. The limiting process obtained by letting a go to zero but keeping 



18 



cr > fjc is no longer a deterministic function, but a "Bernoulli" process, choosing between 
the left and the right potential well with probability 1/2 at integer times (Figure ^). 

Another difference lies in the distribution of barrier crossing times in the transition 
window. In the asymmetric case, paths may overcome the saddle as soon as t ^ — cj^/^, 
and are unlikely to return to the shallower well. In the symmetric case, paths may jump 
back and forth between both wells up to time ^/a, before settling for a potential well. 



3.4 Modulated noise intensity 

Other mechanisms leading to stochastic resonance have been examined, for instance peri- 
odic forcing which is not deterministic, but affects the noise intensity, a situation arising 



in power amplifiers |D&j| . This case can be analyzed by the same method as the previous 
ones, but the results are quite different. 

Let us consider the motion in a static symmetric double-well potential, described by 
the SDE 

dxt = - [fiQXt - X?] dt + ^g{t) dWu (3.14) 

£ ye 

where > is fixed and g{t) is periodic. We assume that g{t) ^ e|logcr|//io foi" all 



t. Theorem shows that for sufficiently weak noise, paths starting at time to at the 
bottom ^/JIQ of the right-hand potential well remain concentrated in a strip around yT^i 
with width proportional to ag{t) / This holds as long as the spreading is smaller 
than a constant times the distance yjjio between well and saddle. The probability to cross 
the saddle before time t is bounded by 

P(t) = C7(t,e)exp|-^f|^^ |, where g{t) = sup g{s) (3.15) 



and K > 0. Thus \i g reaches its maximum ^(0) at time 0, the probability to see a transition 
during one period satisfies 

Ptrans ^ C(l, s) e-'^'^^/^^^ for ^ ^ = (3.16) 

Note that here, as the potential is static, the threshold value for the noise intensity can 
be guessed from Kramers' time, assuming constant g. Taking into account that g is not 
necessarily constant, we see that for a > cr^, transitions are likely to happen in the time 
interval during which ag{t) ^ 2^o- Tor instance, if g{t) behaves quadratically near its 
unique maximum, this transition window is given by 

< const 5(0) [l-'^y (3.17) 

In contrast to the previous cases, however, the transition times are less concentrated, in 
the sense that for times to < < ^2 before the transition window, 

P{ti) ~ P(t2)®(*2)/9(ti))'. (3.18) 
A similar argument as in the previous cases shows that for o" > (Tc, 

1 f 2/ioA 1 

^trans ^ o " ^onst exp<^ -K— r ), (3.19) 

2 1^ e|logo-|J 
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where A is the length of the transition window. 

If the potential is made asymmetric, so that a constant term Aq > is added to the 
drift term in (3.14), the critical noise intensities needed to reach the saddle from the 
shallower left-hand well and the deeper right-hand well will be different (one can check 
that their ratio is 1 + 0{\q/ )). As a consequence, transitions from the shallower to 
the deeper well will be likely as soon as the noise intensity exceeds the smaller threshold, 
while transitions in both directions are likely when the noise intensity exceeds the larger 
threshold. When the noise intensity drops below the smaller threshold again, xt will be 
in the deeper right-hand well with larger probability. The net effect is that xt will visit 
both wells near integer times, but has only small probability to remain in the shallow well 
after the transition window. Thus the periodic signal is not amplified in the same way as 
discussed before, but nevertheless we observe an amplification mechanism which allows to 
read off at which times the threshold is exceeded. 



4 Hysteresis 

Hysteresis is another characteristic phenomenon of bistable systems. Let us consider again 
the motion in the Ginzburg-Landau potential ( |3.1| ) with fi = 1 and X{t) = —Acos(2TTt), 
but without imposing the restriction A < X^. In the deterministic case a = the equation 
of motion reads 

s^ = xt-xl + Xit). (4.1) 

We may ask the question: How does xt behave, as a function of X{t), in the adiabatic 
limit e — > 0? Intuitively, xt will always track the bottom of a potential well. A naive way 
to see this is to set formally e equal to zero in (|4.l| ): We obtain the algebraic equation 
X — x^ + X = 0, which admits three branches of solutions (see Figure]^; X^{X) and X* (A) 
correspond to potential wells, and -^o(A) to a saddle, which exists only for |A| < Ac. 

If the amplitude A is smaller than the critical value Ac, there are always two potential 
wells separated by a barrier. Hence xt will always track the bottom of the same well in 
the limit e — > 0, so that the instantaneous value of A is sufficient to determine the state 
(provided we know in which potential well the process started). 

If A is larger than Ac, however, a saddle-node bifurcation point is crossed whenever 
|A| reaches Ac from below: The potential well tracked by xt disappears, so that xt jumps 
to the other well, which is unaffected by the bifurcation (Figure ||). As a result, the 
state Xt is not uniquely defined by the instantaneous value of A if |A| ^ Ac: xt tracks 
the bottom X* (A(t)) of the left-hand well if A increases, and the bottom X^{X{t)) of the 
right-hand well if A decreases. This phenomenon is called hysteresis. The hysteresis cycle 
consists of the branches X^{X), \X\ ^ A, and two vertical lines on which |A| = Ac. It 
encloses an area = 3/2, called static hysteresis area. In many applications, x and A 
are thermodynamically conjugated variables, and the hysteresis area represents the energy 
dissipation per cycle. 



4.1 Dynamical hysteresis and scaling laws 

Consider now what happens when e is small but positive. The solutions of Equation (4.1) 
will not react instantaneously to changes in the potential, so that the shape of hysteresis 
cycles is modified. It is important to understand the e-dependence of quantities such as 
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Figure 6. Periodic solutions of the deterministic equation (4.1), (a) in a case where the 
ampHtude A of \{t) is smaUer than Ac, and (b) in a case where it is larger than Ac. The 
enclosed area scales like eA in case (a), and like Aq -\-£^^^{A— Ac)^^'^ in case (b), where 
is the static hysteresis area. Potential wells X'^{\) are displayed as full curves, the saddle 
(A) as a broken curve. 



the average of xt over one period, the value of A when xt changes sign, and the area 
enclosed by the hysteresis cycle. It is known that there are constants 71 > 70 > such 
that the fohowing properties hold: 

• If A ^ Ac + 7oe, then xt cannot change sign (except during the very first period, if 
the process does not start near the bottom of a well). There exist two stable periodic 
orbits, one tracking each potential well (Figure ^). Each encloses an area A satisfying 

A X Ae, (4.2) 

and the average of xt over each cycle is nonzero. The notation x is a shorthand to 
indicate that c^Ae ^ A^ c+Ae for some constants c± > independent of e and A. 

• If A ^ Ac + 7ie, then xt changes sign twice per period. All orbits are attracted by the 
same periodic orbit (Figure ^3), corresponding to a hysteresis cycle with zero average 
and area A satisfying 

A-Ao>^e^/^{A-X,)^/\ (4.3) 

When xt changes sign, the parameter A satisfies |A| — Ac x e^/^(A — Xc)^^^. 

• If Ac + 7o£ < A < Ac + 7ie, several hysteresis cycles may coexist, some of them 



satisfying (4^) and others satisfying ([4.3|). 



The scaling law (4.3) was first derived in pGRM| ] for A—Xc of order 1, where Equation ( |4.1| 



was used to model a bistable laser. The case where A is close to Ac has been analysed in 



An equation qualitatively similar to (|4.1| ) describes the dynamics of a Curie -Weiss 



model of a ferromagnet, subject to a periodic magnetic field X{t), in the limit of infinitely 
many spins | Mar| . The transition between the small and large amplitude regimes has been 



called "dynamic phase transition" in [TO]. 



The magnetization obeys a deterministic differential equation only in the limit of infi- 
nite system size. The effect of the number of spins being finite can be modeled, in first 



approximation, by an additive white noise of intensity proportional to 1/vN |Mar]. It is 
thus of major importance to understand the effect of additive noise on the properties of 
hysteresis cycles. 
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Langevin equations have already been studied for multi-dimensional Ginzburg-Landau 
potentials. Then, however, the mechanism leading to hysteresis is different, because there 
is no potential barrier between stable states. Numerical simulations |RKP] suggested that 
the area of hysteresis cycles should follow the scaling law A x while various 

theoretical arguments indicate that A ^ e^l'^A^I'^ [DT, SD, ZZ]. It is not clear whether 
such a scaling law exists for the Ising model |^RN |. 



For clarity, we will keep interpreting xt as magnetization and A (4) as magnetic field. 
There exist, however, many other instances where the dynamics is described by a peri- 
odically forced Langevin equation. For instance, in models for the Atlantic thermohaline 
circulation, xi represents the salinity difference between high and middle latitude, and 
A(i) represents the atmospheric freshwater flux Ra|. The effect of additive noise on 



this system has been investigated, for instance in |Ce|, while the properties of hysteresis 



cycles were considered in particular in [Mc]. 

The fact that additive noise may create relaxation oscillations has been discussed 



in [Fr2|, where the motion of a light particle in a randomly perturbed field is investigated 
with the help of large deviation theory. 



4.2 The effect of additive noise 

We consider the Langevin equation 

1 



dxt = -\xt-x\- A cos(27rt)] dt + dW^, 
e ye 



(4.4) 



where A > 0. We denote ^ — Ac by ao, but in contrast to Section ^ (where ao had opposite 
sign), we do not impose that ao is a small parameter, and we allow positive as well as 
negative oq. Let us fix a deterministic initial condition (to = — l/2,xo > 0), such that 



the solution x^''* of the deterministic equation (4.1) with 
right-hand potential well. 

We are interested in the quantity 



„det 



A{e,o) 



1/2 
1/2 



xtA'(t)dt 



xq is attracted by the 



(4.5) 



measuring the area enclosed by xt in the (A, x)-plane during one period (xj does not 
necessarily form a closed loop, but A still represents the energy dissipation). Aie^ 0) is the 
area enclosed by xf*^*, and behaves like (O) or ( |4.3| ). For positive u, ^(£,0") is a random 



variable, the distribution of which we want to characterize. 

Another random quantity of interest is the value A'^ of the magnetic field when Xf 
changes sign for the first time: 



A° = A(r°), 



T° = inf{t > to: Xt ^ O}. 



(4.6) 



Results from Section |3| already allow us to make some predictions for the case ao < 0. 
For o" ^ (Jc = (|ao| V e)^/"^, xj is unlikely to switch between potential wells, so that A(e^ o) 
will be concentrated near the deterministic value ^(e, 0), which is of order e (Figure |8|a). 
For (T 3> (Tc, Xt is likely to cross the potential barrier at a random time which behaves 
typically like — cr^/^. The corresponding field A'^ behaves like Ac — o"^/^ (Figure ^c). Thus 
additive noise of sufficient intensity will lead to a hysteresis area which is smaller, by an 
amount of order o"^/^, than the static area ^o- 
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Figure 7. Definition of the parameter regimes for hysteresis cycles, shown in the plane 
(ao = A — Ac, cr) for a fixed value of e. The behaviour of the hysteresis area A{e, a) in each 
regime is described in Theorem 4.1. Typical hysteresis cycles are shown in Figure @. 



The same behaviour can be shown to hold for positive ao up to order e. In this case, 
the potential barrier vanishes during a short time interval, which is too short, however, 
for xt to notice. For ao > e, there is a similar transition between a small-noise regime 
(Figure §3), where xt is likely to track the deterministic solution, and a large-noise regime, 
where it typically crosses the potential barrier some time before the barrier vanishes. The 
threshold value of a delimiting both regimes is again deduced from the variance of the 
equation linearized around x*^'^*, and turns out to be Uc = {e^/aoY^'^ . For a > ctc, the 
typical value of the field when xt changes sign is again found to behave like Ac — cr^/^ . 

We thus obtain the existence of three distinct parameter regimes, with qualitatively 
different behaviour of typical hysteresis cycles. We summarize the main results in the 
following theorem, and give some additional details afterwards. Many estimates contain 
logarithmic dependencies on ao, cr and e. In order not to overburden notations, we will 
assume that a and ao behave like a power of e (ao may also be a constant), and denote 
I log e I by is- The regimes are those indicated in Figure |7|, but some results are only valid 
if we exclude a logarithmic layer near the boundary, for instance Case II corresponds to 
oo > 71^ and a ^ const (ey^)^/^/^^. 

Theorem 4.1 ([^G|, Theorems 2.3-2.5]). 

• Case I: (Small-amplitude regime) 

The distribution of the random area A{£, a) is concentrated near the deterministic 
value A{£,0) X Ae. There are two subcases to consider: 

— In Case la, A{e, a) can be written as the sum of a Gaussian random variable with 
variance of order a'^s, centred atA{s,0), and a random remainder. The remainder 
has expectation and standard deviation of order a'^ig at most. 

— In Case lb, the distribution of A{£,a) is more spread out. Expectation and stan- 
dard deviation of A{e,a) — A{e,0) are at most of order a'^ie, which may exceed 
A{e,0). 

• Case II: (Large- amplitude regime) 

The distribution of A{£, a) is concentrated near the deterministic value A{£, 0) which 
satisfies (|4.3[) . 
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Figure 8. Typical random hysteresis "cycles" in the three parameter regimes of Figure 0. 
Deterministic solutions are shown for comparison, (a) Case I, small-amplitude regime 
(here e — 0.05, oo ~ —0.1, cr — 0.025): Paths typically stay close to the deterministic 
cycle, which tracks a potential well, (b) Case II, large-amplitude regime (here e = 0.02, 
ao = 0.1, a = 0.05): Typical paths are close to the deterministic cycle, which switches 
between potential wells, (c) Case III, large-noise regime (here e = 0.001, oq = 0, tr = 0.16): 
Paths are likely to cross the potential barrier when |A| is of order Ac — a^^^ . 



— In Case Ila, A{e, a) can he written as the sum of a Gaussian random variable 
with variance of order (T^(ey^)^/'^, centred at A{£,0), and a random remainder. 
The remainder has expectation and standard deviation of order a'^£^{e^/ao)~'^^^ at 
most. 

— In Case lib, we can only show that the distribution of A{£, a) is concentrated in 
an interval of width (ey^)^/'^^^ around A{£ , 0) . 

• Case III: (Large-noise regime) 

The distribution of A{£,a) is concentrated near a (deterministic) reference area A 
satisfying A — Aq s< —a^^^. The standard deviation of A{£,(t) is at most of order 
^4/3^2/3^ anrf its expectation belongs to an interval 

[A - 0{a^l^ej^),A + 0(a2£2) + 0{ei,)\ . (4.7) 

In the case where qq ^ e and a ^ Oq''^; ihe term 0{££s) has to be replaced by 
Oiey^lie/a'^/^). In both cases, the distribution decays faster to the right of A than 
to the left. 

In Regimes I and II, the main effect of additive noise is to broaden the distribution of 
the area, which remains concentrated, however, around the corresponding deterministic 
value. In Regime III, on the other hand, the hysteresis area obeys a completely new scaling 
law, which is determined by the noise intensity rather than by frequency and amplitude 
of the driving field. 

The Gaussian behaviour of A in Cases la and Ila is obtained in the following way. The 
deviation yt = xt — xf*^^ from the deterministic solution satisfies an equation of the form 
(2.24) with g = 1, whose solution obeys the integral equation 



yt = ^ [ e"(*'^)/^ dWs + - [ e"(*'^)/^ 6(y„ s) ds, (4.8) 
Jto ^ Jto 
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where b{y, s) = 0{y'^). For small values of yt, the first term dominates the second one. Its 
contribution to A{£, a) — A{e, 0) = — / yu^'{u) du can be written as 



^ 7(to + l,s)dW„ where j{t,s):=-[ e"("'')/^ A'(u) du. (4. 

V ^ J to Js 



<-tQ + l ft 
'to 

The variance of this term is given by 



1 /■* 

a^eTito + l,to), where T{t,to):=— / -f{t,sfds. (4.10) 

^ Jto 

The integral r(to + 1, ^o) depends only on properties of the deterministic solution xf*^^ via 
the curvature a{t). The auxiliary function 7(t,s) can be evaluated by partial integration, 
its leading term behaving like —e\'{s)/\a{s)\. 

In Case I, r(to + l,io) is of order 1, and thus the contribution of the linear term to 
the variance of the area is of order a'^e. In Case la, one can show that the Gaussian term 
dominates the distribution of A{£,cr) near A{e,0), in the sense that 

F{\A{e,a)-A{e,0)\ ^ H} ^e-^'^Vi-'^) (4.11) 

holds for some constants C, k > 0, and for all H smaller than a constant times y^(|ao| V 
e)^/3 if |ao| ^ e^^'^jtj^, and all H smaller than e/4 if |ao| ^ e^l'^jtj'^. Note that the 
upper bound ( [4.11 ) is exponentially small for the maximal value of H ^ except on the upper 



boundary of Region la. 

In Case lb, the Gaussian term no longer dominates, but one can still show that 

P{|^(e,a) -^(e,0)| ^ U] ^ ^ g-^-f^/^'^'^-) (4.12) 



up to = const |aop^^4. Again, (4.12) is exponentially small except on the upper 
boundary of Region lb. 



Estimates ( [4.11| ) and ( 4.12| ) control the tails of the distribution of A(e^a^ in a neigh- 



bourhood of ^(e,0). The quartic growth of the potential Vix^t) for large implies, on 
the other hand, that 

F{\A{e,a)-A{eM ^ H} ^'^e'^^''^^ (4.13) 

for all H larger than some constant (of order 1). In fact, this estimate holds in all 
parameter regimes, since it does not depend on the details of the potential near x = 0. 
This still leaves a gap between the domains of validity of ( [4.11| ) and ( 4.12| ), and of ( [4.13 ), 



which is due to the existence of a second potential well. In fact, the distribution of the 
hysteresis area will not be unimodal. Sample paths are unlikely to cross the potential 
barrier, but if they do so, then most probably near the instants of minimal barrier height, 
in which case they enclose an area of order 1. Hence the density of ^(£,0") will have 
a large peak near ^(e, 0), and a small peak near areas of order 1 (more precisely, near 
!^j,{XX{X) - X^_{X)) dA), see Figure |. 

In Case Ila, the distribution of ^(e, cr) near ^(e, 0) is again dominated by a Gaussian, 
stemming from the linearization of the SDE around j;^*^*. The integral r(to + l,to) ™- 
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Figure 9. Sketches of the distribution of the hysteresis area A{e, a) in the different param- 
eter regimes. Regime la: The area is concentrated near the deterministic area A{£, 0) x Ae. 
There is a small probability to observe areas of order 1. Regime lb: The distribution is 
more spread out. Regime Ila: The area is concentrated near the deterministic area A{e, 0) 
of order ^0 + (£\/ao)^^'^- Regime lib: The area is concentrated in an interval of width 
(ey/oo)^^^ around Ao + {e^/oo)'^^^ , but we do not control the distribution in this interval. 
The broken curve shows an extrapolation of the estimates outside this interval. Regime III: 
The area is concentrated near A, which is of order Ao — a^l"^ . The distribution decays 
faster to the right than to the left. 



,2/3 



( 4.10 ) is found to behave like e ^/^Oq''^, leading to a variance of order (T^(ey^)^/^ and to 



the bound 



P{|^(e,(7) -Aie^ 



>H}^ -e 



(4.14) 



valid for H smaller than a constant times £^/ao. 

Unfortunately, this estimate cannot be extended to Case lib. The reason is that during 
the jump of xf'^^ to the left-hand potential well, a zone of instability is crossed where paths 
are strongly dispersed. The maximal spreading of paths is of order (T(e-y/ao)~^/®, which 
is too large, in this regime, to allow for a precise control of the effect of nonlinear terms. 
This does not exclude the possibility that the bound ( 4.14| ) remains valid for larger a. 

However, the value A*^ of the magnetic field when xt changes sign can be described in 
all of Regime II. One can show that 

P{|A0| < Ac-i} ^ |e-'^(l^l'''^^v^)/'^' for -Li(e^)2/3 ^ L ^ Lo/4 (4.15) 

P{|A0| > Ac + L} ^ 3g-.L/(a2(ev/^)3/24) 



for 

for L ^ L2{£^f/^. 



(4.16) 
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Note that [A*^! cannot exceed Ac + ao- These bounds mean that the distribution of A*^ is 
concentrated around the deterministic value of order Ac + (ey^)^/^, and decays faster to 
the right than to the left. They can be used to show that in Case lib, 

P{^(e, a) - A{£, 0)^-H}^^ ^-kH^i^Io^ ^ (417) 

P{^(e, a) - A[e, 0) ^ +i7 } ^ ^ ^-K{e^,Yi^HI{aH,) ^^ -^^^ 

for (e^)2/34 ^ i7 ^ {e^fa^Yl'^iz. Thus the probability that ^(e, cr) - A{e, 0) is outside 
an interval of size (ey^)^/^^^ is very small. We do not control, however, what happens 
inside this interval. 

In Case III, the large-noise regime, most sample paths are driven over the potential 
barrier as soon as the magnetic field reaches a value of order Ac — cr^/^. Rather than 
comparing Aie^ a) to its deterministic value, we should compare it to a reference area A 
given by 



U= r xf''+{-X'{s))ds+ ['^\f''-{-X'{s))ds, (4.19) 

^ 7-1/4 Jti 



where xf*^*'^ are solutions tracking respectively the right and left potential well, and ti x 
-^2/3 is the typical jump time. Checking that A — Ao scales like -^4/3 is straightforward . 

The probability of deviations of A{e, a) from A can be estimated by bounding sep- 
arately the integrals between —1/4 and ti and between ti and 1/4. The results differ 
slightly in two regimes. If oq ^ e or cj > Oq^^, then 

P{^(e, a)-A^-H}^- ^-^"'1^' +^ (4.20) 

P{^(e, a)-A^ +H} ^ - e-^^/('^'^-) +- g-'^^/^^^-) (4.21) 

holds for some C, k > and all up to a constant times a'^^^le. If oq ^ e and a ^ clq^^, 
two exponents are modified: 

F{A{e, a) - i ^ -H} ^ - e"'^^'^'/'^' +- e-«'^'/(^^^-) (4.22) 
P{^(e, a)-A^+H}i^- e-^/(-'^^) +- ''"^/(^v^^^^) . (4.23) 

Again, the distribution decays faster to the right of A than to the left, guaranteeing that 
^(£,0") is likely to be smaller than the static hysteresis area ^o- The second term on the 
right-hand sides of ( [4.20| ) and ( 4.22 ) does not depend on H: It bounds the probability 



that the paths do not cross the potential barrier, and enclose an area close to zero. The 
situation is opposite to Case I: The distribution of the hysteresis area has a large peak 
near A and a small peak near 0. 

The qualitative behaviour of the distribution of A{e, a) is sketched in Figure When 
crossing the boundary between Regime I and Regime III, the probability of paths crossing 
the potential barrier increases, so that the peak near A{e, 0) shrinks while the peak near ^0 
grows. When approaching the transition line between Regimes III and II, the distribution 
of the area ^(e, cr) becomes more spread out and more symmetric, before concentrating 
again around the large- amplitude deterministic value of ^0 + 
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5 Bifurcation delay 



In the previous section, we had to deal in particular with the slow passage through a 
saddle-node bifurcation. This section is devoted to the slow passage through a (symmetric) 
pitchfork bifurcation. 

We consider again the Ginzburg-Landau potential (|3.1|), but this time with A = 0, and 
a parameter increasing monotonously through zero. As /x changes from negative to 
positive, the potential transforms from a single-well to a double-well potential, a scenario 
known as spontaneous symmetry breaking. In fact, the symmetry of the potential is not 
broken, but the symmetry of the state may be. Solutions tracking initially the potential 
well at X = will choose between one of the new potential wells, but which one of the 
wells is chosen, and at what time, depends strongly on the noise present in the system. 



5.1 Dynamic pitchfork bifurcation 

In the deterministic case (7 = 0, the equation of motion reads 

= fiit)xt - xl (5.1) 
Its solution xf'^^ with initial condition xf'^^ = xq > can be written in the form 



^det 



c(xo,t)e"(*'*o)/% a{t,to)= fi{s)ds, (5.2 



to 



where the function c{xQ,t) is found by substitution into (5.1). Its exact expression is of 



no importance here, it is sufficient to know that < c{xo,t) ^ xq for all t. 

Assume that /i(t) is negative for t < and positive for t > 0. If we start at a time 
to < 0, the solution ( |5.2| ) will be attracted exponentially fast by the stable origin. The 
function a(t, tg) is negative and decreasing for tQ < t < 0, which implies in particular that 
Xq''* is exponentially small. For t > 0, the function a{t,to) is increasing, but it remains 
negative for some time. As a consequence, Xj**^* remains close to the saddle up to the first 
time t = n(to) for which a{t, to) reaches again (if such a time exists). Shortly after time 
n(to), the solution will jump to the potential well at fJ-it), unless xq is exponentially 
small. n(to) is called bifurcation delay, and depends only on fi and to. For instance, if 
fi{t) = t, then a{t,to) = ^{t^ - tg) and n(to) = |to|- 

The existence of a bifurcation delay may have undesired consequences. Assume for 
instance that we want to determine the bifurcation diagram of x = /ix — x^ experimentally. 
Instead of measuring the asymptotic value of xt for many different values of /i, which is 
time-consuming (especially near = where xj decays only like l/\/i), one may be 
tempted to vary slowly during the experiment. This, however, will fail to reveal part of 
the stable equilibrium branches because of the bifurcation delay. A similar phenomenon 
exists for the Hopf bifurcation | Nelj , Ne2| ] . 



The delay is due to the fact that xt approaches the origin exponentially closely. Noise of 
sufficient intensity will help driving the particle away from the saddle, and should therefore 
reduce the bifurcation delay. The obvious question is thus: How does the delay depend 
on the noise intensity? 

In the case of a quadratic potential V{x,t) = — i/i(t)x^, this question has been in- 
vestigated and compared with experiments, with the result that the delay behaves like 
V^logoi |SMC| , |SHA| . Similar results were obtained in [|JL|] . Our techniques allow 



us to derive rigorous bounds for the nonlinear equation. 
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5.2 Effect of additive noise 

We consider the nonlinear SDE 

dxt = -[fi{t)xt-xf]dt + ^dWt. (5.3) 

£ ye 

The results presented here are a particular case of those obtained in [ |BG1 ] , which apply to 
more general nonlinearities. We assume that = t + 0{t'^) is monotonously increasing 
on an interval [to,T] or [to,C)o), to < 0. We denote by xf*^^ and xt the solutions of the 
deterministic and stochastic equations with given initial condition xq > 0. 

We already know that xf*^^ decreases exponentially fast for t < 0. Results from Sec- 
tion § show that paths are concentrated in a neighbourhood of order a of xf^^ on any 
time interval [ioj^i] bounded away from zero, so that we only need to worry about what 
happens after time ti, when x'^'^^ is already exponentially small. 

The dispersion of paths will be controlled by the variance-like function 

v{t) = vo e2"(*'*^)/- +- /* e2"(*'^)/- ds, (5.4) 



where vq is a positive constant. One can show that this function grows like a'^/\fj,{t)\ for 
t ^ and remains of order j^/e up to time \/e. Only after time \/e, v(i) grows 

exponentially fast. In analogy with ( p. 16 ), we define a strip 



Bih) = {(x,t): ^ t ^ ^/e, |x -xf*| < h^f^}. (5.5) 

In order to describe the behaviour for t ^ y^, we further introduce the domain 

V{q) = {{x,t): V^^t^T,\x\ ^ y'il- Q)fi{t)}, (5.6) 

where ^> is a parameter in [0,2/3). Note that P(2/3) contains those points in space-time 
where the potential is concave, while 25(0) contains the points located between the bottoms 
of the wells. 



Theorem 5.1 ( [PG1| , Theorems 2.10-2.12]). 

• There is a constant ho > such that for all h ^ hQ^/e/a, the first-exit time Tj2(h) of 
xt from B{h) satisfies 

P*^'"*i{rB(h)< V^I^S^e-'''^', (5.7) 

where 

C. = MVM4±^ and .= '--0{V-e)-Or^). (5.8) 

£^ 2 \ e J 

• Assume that a\\oga\^/'^ = 0{y/£). Then for any q G (0,2/3), the first-exit time tx)[q) 
of Xt from V{q) satisfies 



P^'^^{r^,(,) ^t]^ C(t,6) l°g^l 1. , (5.9) 



where 



C{t,e) = const fi{t)(l + ^i^l^V 



(5.10) 
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Figure 10. A sample path of the SDE with fi{t) = t, for e = 0.01 and a = 0.015. A 
deterministic solution is shown for comparison. Up to time y^, the path remains in the set 
B{h) centred at xf'^^, shown here for h = 3. It then leaves the set ^{g) (here g = 2/3) after 
a time of order ^ £|log cr|, after which it remains in a neighbourhood of the deterministic 
solution starting at the same time on the boundary of 'D{g). 



Assume xt leaves V{q) (with 1/2 < < 2/3/ through its upper (lower) boundary. Let 
x^*^*'^ he the deterministic solution starting at time r = Tx)(g) on the upper (lower) 
boundary of D{q). Then x^'^*'^ approaches the equilibrium branch at \J fi{t) (resp., 
— ■\/li{t) ) like £/ii{t)'^/'^ + \/ ^i{t) e~''"(*''^)/'^, where rj = 2 — 2>q. Moreover, xt is likely 
to stay in a strip centred at x^^*'^, with width of order a j \J [i(t), at least up to times 
of order 1. 



The bound (^.T]), which is proved in a similar way as Theorem 2.2, shows that paths are 
unlikely to leave the strip BiK) if 1 ^ /i ^ hi^^ja. If a is smaller than y^, paths remain 
concentrated in a neighbourhood of the origin up to time -y/e, with a typical spreading 
growing like o j \J\li{t)\ for t ^ — \/e, and remaining of order a/e^^^ for |t| ^ y/e. This is 
again a dynamical effect: Although there is a saddle at the origin for positive times, its 
curvature is so small that paths do not have time to escape before t = y^, see Figure [To| . 

Relation ( |5.9[ ) yields an upper bound on the typical time needed to leave I'(^>), and 
thus enter a region where the potential is convex. Since a(t, y/e) grows like for small 
t, the probability not to leave T>{q) before time t becomes small as soon as 



t > y^e|logcj|. (5.11) 

The last part of the theorem implies that another time span of the same order is needed 
for paths to concentrate again, around an adiabatic solution tracking the bottom of the 
well (at a distance of order e///(t)^/^). One can thus say that the typical bifurcation delay 
time of the dynamical pitchfork bifurcation with noise is of order y^e|log a\. 
As a consequence, we can distinguish three parameter regimes (Figure pi]): 
I. Exponentially small noise: a ^ e-^l^ for some K > Q. 
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Figure 11. Depending on the value of noise intensity a and drift velocity e, the random 
bifurcation delay has one of three qualitatively different behaviours. If a is exponentially 
small in 1/e, the typical delay is macroscopic. For a larger than y^, there is no such delay, 
and paths make large excursions near the bifurcation point. In the intermediate regime, 
the typical delay is of order ^e|logcr|. 



At time -y/e, the spreading of paths is still exponentially small. In fact, one can extend 
Relation ( ^ ) to all times for which a{t,^/e) < K. li K is larger than a(n(to),0), 
where n(to) is the deterministic bifurcation delay, most paths will track the determin- 
istic solution and follow it into the right-hand potential well. 
II. Moderate noise: ^ u <C ^/e for some p > 1. 

The bifurcation delay lies between ^/£ and a constant times y^ejlog a\ ^ ^(i^p)/^ with 
high probability. One can thus speak of a "microscopic" bifurcation delay. 
III. Large noise: a ^ y/e. 

The spreading of paths grows like a / ■\/|/u(t)| at least up to time —a. As t approaches 
the bifurcation time 0, the bottom of the potential well becomes so flat that the 
paths are no longer localized near the origin and may switch wells several times before 
eventually settling for a well. So for large noise intensities, the concept of bifurcation 
delay should be replaced by two variables, namely the first-exit time from a suitably 
chosen neighbourhood of the saddle and the time when the potential wells become 
attractive enough to counteract the diffusion. 

One can also estimate the probability to reach the right-hand potential well rather 
than the left-hand potential well. Loosely speaking, if Xs reaches before time t, it has 
equal probability to choose either potential well. It follows that 

P*0'^o{xj ^ 0} = ^ + ip*0'^o{xs > Vs G [tQ,t\]. (5.12) 

One can show that for t = the second term on the right-hand side is of order 

^^e-l"(0■*°)l/^ (5.13) 
a 

and for larger t, this term will be even smaller. Thus in case II, paths will choose one 
potential well or the other with a probability exponentially close to 1/2. 

In the particular case where /i(t) > const t for all t > 0, the height of the potential 
barrier grows without bound. This implies that once xt has chosen a potential well, its 
probability ever to cross the saddle again is of order e"'^""** . 

The existence of three parameter regimes has some interesting consequences on the 
experimental determination of a bifurcation diagram. Assume we want to determine the 
stable equilibrium branches by sweeping the parameter with speed e. Regime II is the 
most favourable: In Regime I, part of the stable branches cannot be seen due to the 
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bifurcation delay, while in Regime III, noise will blur the bifurcation diagram. For a given 
noise intensity cr, the sweeping rate e should thus satisfy 

^2 «£« (l/|loga|)'/^ (5.14) 

in order to produce a good image of the stable equilibria. As long as o"^ <^ l/\loga\^^^ , 
increasing artificially the noise level allows to work with higher sweeping rates, but of 
course the image will be more and more blurred. 

On the other hand, the relation between noise and delay can be used to measure the 
intensity of noise present in the system. If a bifurcation delay is observed for a sweeping 
rate Eq, repeating the experiment with slower and slower sweeping rates should ultimately 
suppress the delay. If this happens for e = ei, then the noise intensity is of order e~™"^* Z^^. 



6 Generalizations 

6.1 Multidimensional systems 

We now return to n-dimensional equations such as ( |1.11| ). As in Section |2|, we start by 
examining the linear equation 

dyt = -A{t)ytdt + ^Git)dWt, yo = 0, (6.1) 

obtained, for instance, by linearizing the equation around a given deterministic solution 
(in that case, the matrices A and G may depend on e). If the drift term derives from a 
potential, then A is necessarily symmetric, but we will not impose such a restriction here. 
We will assume, however, that we are in the stable case, that is, all eigenvalues of j4(t) 
have real parts smaller than some constant — oq < 0. 
The solution of ( |6.1| ) can be written as 

yt = ^ f Uit,s)G{s)dWs, (6.2) 
Jo 

where U{t,s) = U{t)U{s)-^ is the propagator of the deterministic equation ey = A{t)y, 
and we denote by U{t) its principal solution, i.e., £U{t) = A{t)U(t) and U{0) = 1. The 
random variable yt has a Gaussian distribution, with zero expectation and covariance 
matrix 

Cov{yt} = — U{t,s)G{s)G{sfU{t,sfds. (6.3) 
£ Jo 

In order to determine the asymptotic behaviour of the covariance matrix, we first note 
that it is a particular solution of the ODE 

£-^X{t) = A{t)X{t) + X{t)A{tf + a'^G{t)G{tf. (6.4) 



The general solution of ( |6.4D can be written as 

X{t) = X{t) + Uit) [X{0) - X{0)] U{tf, (6.5) 
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where X{t) is a particular solution admitting an asymptotic expansion (c. f. [ |Wa| ]) 

X{t) = Xoit) + eXiit) + e^X2{t) + ■■■ (6.6) 
The terms of this expansion satisfy the Liapunov equations 

A{t)Xo{t) + Xo{t)A{tf = -a^G{t)G{tf (6.7) 

A{t)Xn{t) + Xn{t)A{tf = ^^Xn-i{t) Vn ^ 1. (6.8) 

at 



The solution of (|6.7] ) admits the integral representation (see for instance [Be]]) 

Xoit) = / e^(*)^ G(t)G(t)^ e^^^^ ds, (6.9) 



and similarly for the solutions of (|6.^). Note that Xo corresponds to the asymptotic 
covariance matrix of ( |6.2[) if A and G are constant in time 

At any fixed time t, the distribution of yt is concentrated in an ellipsoid of the form 
{y , Cov {yt}~^ y) ^ const, provided the (Gaussian) distribution of yt is nondegenerate. 
Assume for the moment that G(t)G(t)"^ is uniformly positive definite. Then Cov{yt}~^ 
exists and so does X(s)~^. In addition, we find that (T^X(s)~^ is bounded in norm. The 
following result, which generalizes Proposition from the one-dimensional case, shows 
that the paths {yt}t^o are concentrated in sets {(t,y): (y, X(t)~^y) ^ const}. The proof 
will be given in pG4| ]. 

Theorem 6.1. Assume that a'^Xis)~^ is bounded in norm. Then for all t > 0, h > 
and any k £ (0, 1/2), 



!)0,0 



sup {y,,Xis)~'ys) ^ h^} < C„(t, e) e-'^'^'^^-^^^", (6.10) 



where 



a(t,.) = (^ + l)(^) . (6.11) 



The exponential growth of the prefactor as a function of the dimension implies that 
the estimate ( |6.10D is only useful for h > yjn. This dependence, however, is to be expected 
as the tails of n-dimensional (standard) Gaussians only show their typical decay outside 
a ball of radius proportional to yjn. 

Under the condition that cj~^X(s) and (T^X(s)~^ are bounded in norm, one can show 
that a bound similar to ( |6.10 ) holds if yg obeys a nonlinear perturbation of Equation (6.1), 
for all h smaller than a constant times . 



6.2 Coloured noise 

White noise has no time-correlations, and is thus appropriate to model the random influ- 
ence of a "fast" environment on a "slow" system, if the relaxation time of the fast system 
is negligible. If the relaxation time is short but not negligible, one has to use coloured 
noise instead of white noise. 

The simplest (still Markovian) model of coloured noise is the Ornstein-Uhlenbeck 
process 

Zt = a [\-^^'-''>dWs, (6.12) 
Jo 
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whose autocorrelation function ^.{ZsZt} = (cr^/27) e"'''^*"'*) [1 — e"^'''*], s < t, decays ex- 
ponentially in t — s. 

Let us examine the influence of such coloured noise on a slowly time-dependent, one- 
dimensional, linear system, described by the SDE 

dxt = a{et)xt dt + g{et) dZf. (6.13) 

We assume that a{£t) is negative, g{£t) is positive, and that both functions are uniformly 
bounded away from zero. Using the differential representation dZt = —^Zt dt + crdVFf of 
the Ornstein-Uhlenbeck process (|6.12| ), and scaling time by a factor e, we obtain that the 
random variable yt = (xt, Zt) obeys a two-dimensional linear SDE of the form (|6.1|), with 



^W=(t^ -7.(tA G{t) = (^f). (6.14) 



-7 / 

For 7 — > 0, we recover the white-noise case, while a — > corresponds to the deterministic 
limit. 



Theorem 3J- shows that paths are concentrated in a tube {ys, X[s)^^ys) ^ const, with 
covariance matrix X{s) admitting the asymptotic series ( |6.6| ). The leading term Xq{s) is 
found from (|6.9D to be 



Xois) = a 



( g{s) \ 

(6.15) 



2(7 + |a(s)|) 2(7 + |a(s)|) 
9{s) 1 



V2(7 + |a(s)|) 27 / 



In particular, the paths {xsjs^o are concentrated in a strip of width ag{s)/ y2{^ + |a(s)|). 
The effect of coloured noise of the form ( |6.12| ) is thus to narrow the distribution of the 
paths, in the same way as if the curvature were increased by an amount 7. 

This result allows us to make some predictions on the effect of noise colour on stochastic 
resonance. We have seen in Section |3.2| that in the standard case of the asymmetric 
double-well potential with additive periodic forcing and white noise, paths switch between 
potential wells when a exceeds the threshold cxc = (ao V e)^/^. This value was obtained by 
comparing the typical spreading of paths to the minimal distance between the deterministic 
solution and the saddle. Repeating this argument in the case of coloured noise given 
by ( |6.12D , we thus expect the threshold noise intensity to be determined by 

a2 = (aoVe)(7V(aoVe)i/2). (6.16) 

This shows in particular that shorter correlation times of the noise require larger noise 
intensities to enable transitions between the wells. 

A similar stabilizing effect of noise colour can be expected when sweeping the bifurca- 
tion parameter through a pitchfork bifurcation as discussed in Section The bifurcation 
delay will lie between the macroscopic delay of the deterministic case and the microscopic 
delay of the white-noise case. The shorter the correlation time of the noise, the larger will 
the delay be. 
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